Diffusion and ballistic transport in clean one-dimensional conductors 
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Diffusion and ballistic transport are normally thought to be mutually exclusive phenomena. The 
reason is that the notion of diffusion entails inelastic scattering processes which typically lead to 
current relaxation. While diffusive transport is the rule in generic many-body systems, ballistic 
transport can be protected in integrable one-dimensional (ID) systems by nontrivial conservation 
laws. The large diffusive response seen experimentally in nearly ideal realizations of the 5=1/2 
ID Heisenberg model is therefore puzzling. Here, we show that diffusion is, in fact, universally 
present in interacting ID systems subject to a periodic lattice potential and question the conjecture 
that integrability necessarily leads to ballistic transport. Furthermore, we demonstrate that even at 
high temperatures current decay can be governed by anomalously large time scales, posing a serious 
challenge for numerical studies. 



I. INTRODUCTION 

For a generic metal at sufficiently high temperatures, 
electric transport is expected to be scattering limited. In 
d spatial dimensions, the signature of diffusive motion is 
the characteristic long-time decay of the autocorrelation 
function (n r (t)n r (0)) ~ t~ d l 2 . Here, n r represents the 
density of a globally conserved quantity ^2 r n r . In very 
clean metals, however, transport can be a subtle issue 
because constants of motion of the system may prevent 
currents from decaying completely, leading to ideal (bal- 
listic) transport. An important role in our understand- 
ing of strongly correlated electrons is played by integrable 
quantum models. Since these models possess an infinite 
number of conserved quantities, one might expect ballis- 
tic transport to be the rule rather than the exception. 1 
Whether or not diffusive behavior is possible at all in such 
systems is indeed an intensely studied 1 " 13 but still open 
question. The unusual dynamics of systems close to an 
integrable point have been investigated by experiments 
on ultracold gases in ID geometries, 14 and is related to 
the question whether these systems thermalize through 
unitary evolution. 15 

In the thermodynamic limit, ballistic transport can 
be defined from the condition that the current-current 
correlation function {J(t)J(fy), where J is the spatial 
integral of the current density operator and the brack- 
ets denote thermal average, does not decay to zero at 
large times. This happens, for example, in a free elec- 
tron gas, where J is proportional to the momentum op- 
erator and therefore conserved in a translationally in- 
variant system. 13 The dc conductivity is then infinite. 
By contrast, real metals are expected to exhibit diffusive 
transport (not necessarily connected with diffusion in the 
autocorrelation function) , meaning that any induced cur- 
rent decays to zero. The conductivity is finite and limited 
by the strongest relaxation process (Matthiessen's rule). 
A third case is realized in systems where the current op- 
erator itself is not conserved but a conserved quantity 




FIG. 1: In a diffusive channel, the conductivity is limited by 
the dominant of the various scattering processes pictured as 
a serial arrangement of resistors. If part of the current is, 
however, protected by a conservation law, a parallel ballistic 
channel for charge transport is opened. 



Q exists which has finite overlap with J. We can then 
write J = J\\+ Ju with J\\ = ({JQ)/{Q 2 ))Q being the 
part which cannot decay, 9 leading to parallel diffusive 
and ballistic channels as indicated in Fig. 1. This idea 
can be generalized to a set of orthogonal conserved quan- 
tities Q n , (QnQm) = (Qn)^,mj an d leads to Mazur's 
inequality 6 ' 16 

Here, L is the system size and T the temperature. The 
Drude weight D measures the weight of the delta-function 
peak in the real part of the optical conductivity at zero 
frequency, cj'(uj) = 27rD5(uj) + (j re g(uj). In principle, both 
D and a reg (uj = 0) can be nonzero. 2 Weak breaking of 
the conservation laws renders the conductivity finite, but 
in this case the channel with the smallest decay rate sets 
a lower bound for the conductivity. 9 

It is important to note that the rhs of equation (1) 
can vanish even if an infinite number of conserved quan- 
tities exists, as we will discuss in the following example. 
Consider the integrable model of spinless fermions (XXZ 
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model) 



N 

H = jJ2i~2 ( c « c «+i + h ' c ) + 



(2) 

Here N is the number of sites, J the hopping ampli- 
tude, c\ annihilates a fermion at site Z, and A is the 
interaction strength. This model is equivalent to the 
anisotropic spin- 1/2 chain and is exactly solvable by 
Bethe ansatz (BA). 17 At half-filling, (m) = 1/2, the 
excitation spectrum is gapless for |A| < 1 and gapped 
for |A| > 1. The current operator is J = YliJh w ^ n 
ji = —iJ(c\c l+1 — cJ +1 c z )/2 as follows from a discretized 
continuity equation. 

The transport properties of the XXZ model have been 
studied for many years. 1 " 13 At zero temperature, the 
Drude weight can be calculated by BA 18 and is found 
to be finite in the gapless and zero in the gapped regime. 
Mazur's inequality can be used to show that D{T) ^ 
away from half-filling at arbitrary temperatures. 6 Re- 
markably, at half-filling the Mazur bound for the Drude 
weight obtained from all local conserved quantities van- 
ishes identically due to particle-hole symmetry. Since 
this is only a lower bound, it does not imply that D 
itself vanishes. However, one can argue 19 that in the 
gapped regime D should remain zero at finite temper- 
atures. The main open question is whether the Drude 
weight is finite at finite temperatures in the half-filled 
gapless case. Since equation (1) is actually an equality 
if all conserved quantities are included, 20 a nonzero D 
at half-filling requires the existence of a nonlocal con- 
servation law. 9 D(T > 0) 7^ at half-filling has been 
found in two independent BA calculations. 4 ' 10 However, 
these results disagree and they both violate exact rela- 
tions for D(T) at high temperatures. 10 Further evidence 
for D(T > 0) 7^ stems from exact diagonalization 7-9 
and Quantum Monte Carlo (QMC). 5 ' 21 We will discuss 
these numerical works in relation to our own results in 
section III. 

Evidence for diffusion in the spin-spin autocorrelation 
function at high temperatures has been sought via ex- 
act diagonalization, 11 QMC 22 and density matrix renor- 
malization group (DMRG). 12 ' 23 The results at infinite 
temperature seemed consistent with an algebraic decay 
(Sf(t)Sf(O)) ~ t~ a with exponent a close to 1/2 as 
expected for d = 1. At low temperatures, the diffu- 
sive contribution was practically undetectable. 12 Mean- 
while, nuclear magnetic resonance (NMR) 24 and muon 
spin relaxation 25 experiments even found evidence for 
low-temperature diffusive behavior in S =1/2 Heisenberg 
chains, but have so far remained unexplained. 



II. FIELD THEORY AND EXPERIMENT 

In the NMR experiment on the spin chain compound 
Sr2Cu03, spin diffusion is observed as a characteristic 
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FIG. 2: Experimental data for the spin-lattice relaxation rate 
of the spin chain compound Sr2Cu03 at h = 9 T taken 
from Ref. [24] (dots) compared to our theory (blue solid line). 
Without diffusion, 7 = 0, l/(TiT) would be almost constant 
(red dashed line). The theoretical result for 1/Ti depends 
on the experimentally determined parameters J and A. The 
curves shown are obtained for J = 2000 K and A as given in 
appendix D. 



magnetic field dependence of the spin lattice relaxation 
rate, 1/Ti ~ 1/V^. 24 Here, only excitations with momen- 
tum q ~ 0, relevant for the studied transport properties, 
contribute. Clearly, S^CuOa is not exactly an integrable 
system. However, the behavior is exptected to be differ- 
ent depending on whether the diffusion constant is de- 
termined by intrinsic umklapp scattering within the inte- 
grable model or by integrability-breaking perturbations. 
The spin excitations propagating in a given channel only 
contribute to the diffusive response at frequencies which 
are small compared to the relaxation rate in that chan- 
nel. If the Drude weight of the XXZ model is large in the 
regime /i « T « J, then we expect a large fraction of the 
excitations in S^CuOs to propagate in a quasi-ballistic 
channel with a very small relaxation rate. The diffusive 
response should therefore be suppressed compared to the 
case where the integrable model has a dominant diffusive 
channel. 

We now calculate 1/Ti by a standard field theory ap- 
proach based on the Luttinger model 17 assuming that 
there is no unknown nonlocal conservation law. For the 
pure Luttinger model, 1/Ti ~ T in the limit T -> 0. 26 A 
finite decay rate 7(T) is obtained at second order in the 
umklapp scattering amplitude as displayed in appendix 
B and appendix D. For the Heisenberg model (A = 1) 
at half-filling, we find 7 - Tj In 2 (J/T). This implies 



spin diffusion in the sense that (ni(t)ni(0)) ~ T^/^/t at 
large times. At high temperatures such that 7 ^> h but 
still T « J, we find 1/Ti ~ T^jjh. A comparison 
of the essentially parameter-free calculated temperature 
dependence with experiment is shown in Fig. 2. The 
good agreement indicates that a large diffusive response 
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is present in the integrable model near half-filling. Fur- 
thermore, this result shows that umklapp scattering is 
a "dangerously irrelevant" perturbation of the Luttinger 
model, 26 completely changing the behavior of 1/Ti in 
the regime h <C T from a constant to a square-root di- 
vergence 1/y/h, as seen in experiment. 

Our field theory calculation assumed D(T > 0) = 0. 
In this case, the optical conductivity is a Lorentzian with 
width set by 7. If conservation laws protecting the Drude 
weight are present, they can be naturally incorporated 
using the memory matrix formalism. 2 But in this case 
spectral weight would be shifted from the Lorentzian into 
a ballistic part that does not contribute to the tempera- 
ture dependence of 1/(T\T) (see equation (F18)). That 
the experimental points in Fig. 2 are actually mostly 
above the theoretical prediction suggests that D is rather 
small near half-filling. Additional sources for diffusion, 
as for example phonons, might explain why our theory 
does not fully account for the temperature dependence 
ofl/CZiT). 

III. NUMERICAL RESULTS 

In order to clarify the contradiction with previous 
studies that supported a large Drude weight at half- 
filling, 7 ' 8 we used a DMRG algorithm 12 ' 23 to calculate 
the current-current correlation directly in the thermo- 
dynamic limit. Fig. 3(a) shows numerical results for 
(ji(t)jo(0)) at T = 00. By summing over sites we 
obtain C(t) = (J(t)J)/L which asymptotically yields D 
according to equation (1). Remarkably, the results in 
Fig. 3(b) show that C(t) is nonmonotonic and does not 
converge to an asymptotic value for times up to Jt = 11. 
This is true within the critical as well as the gapped 
regime. We conclude that a large time scale persists at 
T = 00. As stated in [8], exact diagonalization data 
should not be used to speculate about the Drude weight 



in the thermodynamic limit at low temperatures, where a 
large length scale (mean free path) surely exists. Our re- 
sults show that finite size data for the integrable model at 
half-filling are not reliable even at infinite temperature. 

Furthermore, we have verified that the QMC results 5 
cannot resolve the small decay rate j(T) <C T and cannot 
be used to conclude that D is finite. To further support 
that 7(T) is nonzero for T « J, we show C(t)/2JT at 
T = 0.2 J in Fig. 3(c). At intermediate times J" 1 <C 
t < I/7, we expect that C(t)/2JT decays linearly with 
a slope proportional to 7 (see equation (F21)). A linear 
fit in this regime yields values which are consistent with 
our theory (see Fig. 3(d)). We also note that the values of 
C(t)/2JT for Jt w 6 are already smaller than the Drude 
weight found in [10] by BA. 

IV. CONCLUSIONS 

To summarize, we have shown that in integrable ID 
systems diffusion can coexist with ballistic transport, in 
the sense illustrated in Fig. 1. This is the scenario for 
the XXZ model away from half-filling. For the half-filled 
case, however, we have argued that the large diffusive re- 
sponse measured experimentally in spin chains and seen 
in our numerical calculations suggests that, contrary to 
common belief, the low-temperature Drude weight is ei- 
ther zero or surprisingly small. 
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APPENDIX A: LOW ENERGY EFFECTIVE MODEL 

Bosonization of the XXZ model, equation (2), in the gapless regime at half-filling leads to the effective 
Hamiltonian 17 ' 27 ' 28 

H = H + H u + i^bc, 

#o = \ I dx [n 2 + (d x 4>) 2 ] , 

H u = X J dx cos(v / 8ttK0), (Al) 

H hc = -2irv\ + J dx(d x R ) 2 (d x L ) 2 -27rv\_ J dx [{d x ^ R ) A + (^0 L ) 4 ] . 

Here, Ho is the standard Luttinger model and H u and H^ c are the leading irrelevant perturbations due to umklapp 
scattering and band curvature, respectively. The bosonic field <\> = 4>r + 4>l and its conjugate momentum II obey 
the canonical commutation relation [(fifa),]!^)] = iS(x — x r ). The long- wavelength (q ~ 0) fluctuation part of the 
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FIG. 3: (a) (ji(t)j o (0)) for A = 0.8 and T = oo. The wave fronts travel with a velocity Voo « 1.2J. (b) C{t) = (J(t)J(0))/L 
at infinite temperatures for various A as indicated on the plot. The solid (dashed) lines correspond to A in the critical (gapped) 
regime, respectively, (c) C(t)/2JT at T = 0.2 J for A = 0.6 (blue solid line), A = 0.8 (black solid line), and A = 1.0 (red solid 
line). The dashed lines are linear fits C(t)/2JT — A — BJt for Jt £ [3.5, 7]. (d) The fitted values in comparison to theory. We 
expect A = vK/4tt(1 + b) and B/2A = 7. Note that a small Drude weight is not detectable in this linear regime (see appendix 
F). 



fermion density is related to the bosonic field by rij ~ ^jKj2md x <$>. The spin velocity and the Luttinger parameter K 
are known exactly from Bethe ansatz: 



ttVI - A 2 
2 arccos A 



K 



7r — arccos A 



(A2) 



In this notation, K = 2 at the free fermion point (A = 0) and K = 1 at the isotropic point (A = 1). The amplitudes 
A, A + , and A_ are also known exactly 28 : 



KT{K) sunV/iT) 
7iT(2 - K) 



r i 



2K-2 J 



2^r (l 



K 



2K-2 



2K-2 



A_ = 



1 ttK 
tan 



2tt 2K-2' 

r ( SK \ r 3 ( 1 ^ 

1 L \2K-2 ) L \2K-2 J 



(A3) 



\2ttK 



^ (2X-2) ^ 3 (2K-2) 



5 



APPENDIX B: RETARDED SPIN-SPIN CORRELATION FUNCTION 



In terms of spin-1/2 operators, the XXZ model in a longitudinal field h reads 

N 

h = J2i J (srs? +1 + s i s i+i + AS i s i+i) - »Bhsn, (bi) 

1=1 

where /ib is the Bohr magneton. The spin operators are related to the spinless fermions by the usual Jordan- Wigner 
transformation. In particular, the z component is mapped to the local density of fermions, S? ~ nj — 1/2. We 
are interested in the retarded spin-spin correlation function Xret(q,u), which can be obtained from the Matsubara 
correlation function 

n i it 

x{q^n) = -^Y j e- i ^ [ dre^iSfWSm) (B2) 

by the analytic continuation iuo n — > uo + i0 + . 

In the low-energy limit, we follow ref. [29] and relate the long- wavelength part of the retarded spin-spin correlation 
function to the boson propagator 

Xret(g,<j) _ /^\ret/ , x _ V _ / R oN 

~Kg^~ m {q ' U) -u*-v*q*-Il«*( q ,u,y (B3) 

We calculate the self-energy Il ret (q,uj) = n^ et (g, oj) +n^(g, uj) by perturbation theory to second order in H u and first 
order in H^ c . In the following, we focus on the half-filling case (h = 0). 
The contribution from umklapp scattering reads 

U r u et (q,u;) = AttKvX 2 [F ret (q,u;) - F ret (0,0)] , (B4) 

where 30 

.ret/. A V M T \ „:„/o_r^ r ( U + Vq \ T f UJ - Vq 



with 



sinh 2if (7ru) ~ 7T V ~ 2^' 



where B(x,y) = T(x)T(y)/T(x + y) is the beta function. For K > 1/2, we need a cutoff in the integral I(z) in 
equation (B5). However, the imaginary part of I(z) does not depend on the cutoff scheme used. 30 The expansion of 
equation (B5) for \uj ±vq\ «T yields both a real and an imaginary part for n^ et (^,o;). The calculation of IL™£(q,uj) 
is also standard. By contrast with U 1 ^ t (q,uj), the result for H^(q,u;) is purely real, as band curvature terms do not 
contribute to the decay rate. The end result is 

U ret (q,co) w -2i<ycu - boo 2 + cv 2 q 2 , (B6) 

with parameters 

2 7 = Y 1 T 4K ~ 3 

b = (Y 2 -7 3 )T 4K - 4 + l^ (B7) 

b 2 bi 

c = -(Y 2 + Y 3 )T 4K - 4 z YiT^ 

C 2 Ci 
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B(K,1-2K) 2 , 
^3 = A-^^cot 2 (7r W l/2- ( B8) 



A = 47rXA 2 sin^Trif) — T(l/2 - if)r(if ). 



Here b\ and ci (62 and C2) are the parts stemming from the band curvature (umklapp) terms, respectively. 

At the isotropic point, A = 1, umklapp scattering becomes marginal and logarithmic corrections have to be taken 
into account. In this case we find 

27 = irg 2 T, 

4 32 7T 

Following Lukyanov, 28 the running coupling constant g(T) is determined by the equation 



1 \ng , 

— I = In 

5 2 



2 T 



(BIO) 



where 7 is the Euler constant. We remark that a similar calculation was attempted in ref. [3], but there the imaginary 
part of the self-energy was neglected. 



APPENDIX C: SPIN DIFFUSION 



In this section we consider the long-time behavior of the spin-spin correlation function (Si +x (i)Sf(0)) at half-filling. 
At zero temperature, it is known 31 that the slowest decaying term in the autocorrelation function for < A < 1 is of 
the form 

p -%wt 

(stmm-—-, ( T =o) (ci) 

with W = v and 77 = (K + l)/2. This oscillating term is attributed to q = tt/2 high-energy particle-hole excitations 
with a hole near the bottom of the band and a particle at the Fermi surface, or a particle at the top of the band and 
a hole at the Fermi surface. At T = 0, the low-energy contributions to the spin-spin correlation function decay faster 
than the high-energy contributions. In contrast, numerical studies show that at high temperatures the oscillating 
terms are still present, but the slowest decaying term has pure power-law decay with no oscillations. 11 ' 12 This slow 
decaying term has been interpreted as due to spin diffusion at high temperatures. Here we will show that a diffusive 
term is already present at low temperatures. In the following we use the boson propagator in equation (B3) to 
calculate the q ~ low-energy contribution to the autocorrelation function in the regime T <C J. 
We can write the low-energy, long- wavelength contribution to (Sf +x (t)Si (0)) as 

G(x, t) = -2 | + J g / + J g e<(--) l ™ X _ r flf T \ (C2) 

with Xret(Q^) given by equation (B3). Doing the integral over q first, we find 

G(x, t) = ^-(1 + b)- 1/2 (l + C )" 3/2 / + J d £ exp[-iW- - i(C a + 2i^\x\/v] ^\^^ - (7 -» -7), (C3) 
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where t = t(l + 6) _1/2 , x = x(l + c)" 1 / 2 , 7 = 7(1 + b)' 1 ^ 2 and f = T(l + 6) 1 / 2 . The integral in equation (C3) has 
both pole and branch cut contributions. We assume 7 <C T at low temperatures. In the limit t ^> I/7, we find 

= (1 + 6)" 1/2 (1 + c)- 3 / 2 [G (x,t) + GW (*,*)], (C4) 

where 

j ~ 8tt 2 v 2 I sinh[7rf (t - x/v)} J + 8tt 2 v 2 \ s inh[7rf (t + x/v)} J ( J 

is the Luttinger liquid result in terms of the rescaled variables and 



is the diffusive term. Therefore, for all < A < 1, at finite temperatures and sufficiently long times the autocorrelation 
function becomes dominated by a low- energy term 



(Sf(t)Sf(O)) - 'V 2 ^' (2V0,f»l/ 7 ) (C7) 



with the universal power-law decay t 1 I <1 expected for diffusion in one dimension. 



APPENDIX D: SPIN-LATTICE RELAXATION RATE 



The linear response formula for the spin-lattice relaxation rate is 

k = W §i^)i 25+ "^^)u- ( m ) 

where A(q) is the hyperfine coupling form factor, con is the nuclear magnetic resonance frequency and 

1 r+00 

S+~(q^)\ h = -J2 / dte^iS+^mih (D2) 

is the transverse spin dynamical structure factor. Here = Sf ± zS^ are the raising and lowering spin operators. 
The expression in equation (D2) is to be calculated using Hamiltonian (Bl) in the presence of a magnetic field h. 
We focus on the experimentally relevant Heisenberg point A = 1. We would like to express 1/Ti in terms of the 
longitudinal structure factor 

1 r+00 

SZZ ^) = atE / dte^iSmSHO)) (D3) 

at zero field. The latter is more easily calculated in the field theory since Sf is related to the local density of fermions 
n/, whereas Sj^ have nonlocal representations in terms of Jordan- Wigner fermions. Although the exchange term in 
the Heisenberg model is isotropic, the magnetic field term in equation (Bl) breaks rotational symmetry. As a result, 
5^ (q,w)\h cannot be directly replaced by 2S zz (q, uo)\h at finite field. However, we note that the longitudinal field 
has a trivial effect on Sf (t) 

Sf(t) = e iHt Sfe~ iHt , = e- iuJet e im S±e~ im , (D4) 

where H = H(h = 0) and uo e — fish is the electron magnetic resonance frequency. If we assume in addition that 
T ^> c<; e , the magnetic field dependence in the thermal average can be neglected and we have 

o r+00 

S + -(q,u)\ h / dte^-^{S!{t)Sm)\^ (D5) 

iV , 7/ J -OO 
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where the correlation function is calculated at h — 0. This leads to the expression for the spin-lattice relaxation rate 

k "J ^\ A (Q)\ 2 S zz (q,u; N -u; e )\ . (D6) 



Using uo e ^> con and 



we find in the regime uo e <C T 



S zz (q,-u; e ) 



21m Xret(q^e) 



I _ e u e /T 



1 



2T 



Im Xret(q,Ue). 



For the in-chain oxygen site in Sr 2 Cu03, we have A(q) = Acos(q/2) with 

,, |2 fc B (2C 6 ) 2 + (2C C ) 2 



2/i 



TT 3 k%J 2 



-{91nK) , 



(D7) 



(D8) 



(D9) 



where fc# is the Boltzmann constant, C b,c are the dimensionless components of the hyperfine coupling tensor, 24 
giN^ = 4.74 x 10 -9 eV and J is the exchange coupling measured in Kelvin. 

At low temperatures T <C J, the integral in equation (D8) is dominated by the q ~ mode of Xret(<L <^0- Performing 
the momentum integral using the retarded correlation function in equation (B3) with the appropriate parameters (B9) 
for the isotropic case yields 



k B (2C b y + (2C C Y 



where 



7T^k 2 B J 2 



X X = 1 + 



(ginh) 2 x 




(D10) 



g g 2 , 5g 3 , 3^3 r 



1 + 



8 64 2tt 
^ 3 ^ , V3 



4 32 

In the limit j(T) ^> co e , we obtain the diffusive behavior 



7T 



(Dll) 



1 

jyf 



h(T) 



IT/\n 2 (J/T) 



(D12) 



We note that the temperature dependence of the effective diffusion constant D s = v 2 /2j is different from the one 
proposed in ref. [24]. 

From a fit of the experimental susceptibility data obtained in refs. [24,32] using the full parameter-free field 
theory 28 ' 33 we obtain J = 2000 ± 200 K. The components of the hyperfine coupling tensor, on the other hand, are 
given by 2C b = 95 ± 10 and 2C C = 44 ± 10. 24 To obtain the curve shown in Fig. 3 we used J = 2000 K and 
2C b = 105, and 2C C = 54. The dashed line in Fig. 3 was obtained by setting 7 = in equation (D10) while keeping 
the temperature dependence in the parameters X\ and X^. This would correspond to purely ballistic transport with 
a temperature-dependent Drude weight as suggested by the Bethe ansatz calculation of ref. [10]. The figure shows 
that a large low-temperature Drude weight cannot be consistent with the temperature dependence of the spin-lattice 
relaxation rate. 



APPENDIX E: FINITE MAGNETIC FIELDS AND THE MAZUR BOUND 



The current density operator is defined from the continuity equation for the fermion density 

dm = -i[m,H] = -(ji -ji-i), 



(El) 
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which for the XXZ model yields 

iJ 



Ji = -y( c z c z+i ~ c Ui c i)- ( E2 ) 



Away from half-filling, the integrated current operator J = nas a finite overlap with the local conserved 

quantities of the XXZ model. The simplest nontrivial conserved quantity is 

= *y \( c \-i c i+i - c \+i c i-i) - A ( c l-i c i - 4 c i-i) (ni+i ~ + (ni-i ~ (4 c i+i ~ c \+i c i) - ( E3 ) 

Here Je is the energy current operator as obtained from the continuity equation for the Hamiltonian density. 6 It can 
be verified that Je is conserved in the sense that [Je,H] = 0. According to Mazur's inequality, equation (1), the 
overlap of J with Je provides a lower bound for the Drude weight: 

The Mazur bound in the limit T — > oo was calculated in ref. [6] 

J 4AW(l/4-m 2 ) 

^Mazur = f l + 8 A2(l/4 + m*) (T » J) ' (E5) 

where m = (Sf) = (ni) — 1/2 is the magnetization in the spin model. 

In the gapless phase, the Mazur bound can also be calculated in the low-temperature regime using the field theory 
representations of J and Je- In the continuum limit, the continuity equation becomes 

d t n(x) + d x j(x) = 0. 



Using the bosonized form n(x) = ^ K /2ird x (j), we obtain for effective model (1) 

J r =-\f^ J dxd t 0=-v^ J ^[n-|(A + + 2A_)n 3 -|(-A + + 6A_)n(^(/>) 2 ] . (E6) 

In the following we neglect the correction to the current operator due to band curvature terms and calculate the 
Mazur bound £>Mazur for the Luttinger model Hq: 



J dxU. (E7) 



The energy current operator for the Luttinger model is 

J E ~ -v 2 J dxUd x (j). (E8) 

At half-filling, the overlap vanishes because J and Je have opposite signatures under the particle-hole transformation 
(j) — > — 0, II — > —II. Away from half-filling, a small magnetic field term in equation (Bl) can be absorbed by shifting 
the bosonic field (here we set /ib = 1) 

<t>^<t>+-J^-x. (E9) 



In this case, the conserved quantity becomes 



2 



Je = -v 2 I dxUd x (/)-hv\l — I dxU = J E + hj. (E10) 



2tt 

We calculate the equal time correlations within the Luttinger model and find 

^Mazur - 2TL ( f2) ~ TV^liW 2 (T '^ ^ J) ' (E11) 
w E / i -r 3K y h ) 

We note that in the limit T/h — > the Mazur bound obtained from the overlap with Je saturates the exact zero 
temperature Drude weight D(T = 0) = vK/An. 18 

One can also use the Bethe ansatz to calculate the Mazur bound in equation (E4) exactly. We computed the equal 
time correlations using a numerical solution of the nonlinear equations obtained within the Bethe ansatz formalism 
of refs. [34,35]. In Fig. 4, the numerical Bethe ansatz solution for small values of magnetization m is compared to the 
field theory formula (Ell). Remarkably, the free boson result in equation (Ell) fits well the behavior of the exact 
Mazur bound out to temperatures of order J. 
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FIG. 4: Mazur bound ^Mazur as defined in (E4) calculated exactly using the Bethe ansatz (dots) compared to the field 
theoretical formula (Ell) (lines). 



APPENDIX F: DECAY OF THE CURRENT-CURRENT CORRELATION FUNCTION 



The time-dependent current-current correlation function can be written as 



C(t) = \{J{t)J) 



duo 



2tt 1 



-iut 



o/T 



\m(J- J) re tH, 



(Fl) 



where J7)ret(^) is the retarded current-current correlation function. The latter appears in the Kubo formula for 
the optical conductivity 



,(0,) = - 



{j 2 ) 

LT 



+ (J;J}ret(L0) 



(F2) 



where we used that the expectation value of the kinetic energy operator is given by —{J 2 )/T if the system has 
vanishing superfluid density. In the limit of small q, the Kubo formulated can be rewritten as 3 



Cr(cj) 



K 

2^ 



-iuj 



V). 



If we apply the result from the self-energy approach in equation (B3), we have 

vK iuo 



In this case the real part of the conductivity is 



Rea(uo) = (t'(lj) 



2tt cj 2 - U(u) ' 
27 



vK 



2tt [(1 + 6)cj] 2 + (2 7 ) 2 ' 



(F3) 



(F4) 



(F5) 



The self-energy approach predicts zero Drude weight whenever j(T) is nonzero. This result is inconsistent with 
Mazur 's inequality if there exist conservation laws which protect the Drude weight. It is possible to accommodate 
the existence of nontrivial conservation laws by resorting to the memory matrix formalism of ref. [2]. This approach 
starts from the Kubo formula for a conductivity matrix a which includes not only the current operator J = J\. but 
also "slow modes" J n (n > 2) which have a finite overlap with J: 



LT 



(F6) 



The conductivity in equation (F2) is the cr n component of the conductivity matrix. The idea is that if {J n {t)Jn{fy) 
is a slowly decaying function of time, the projection of J into J n governs the long-time behavior of the current- 
current correlation function and consequently dominates the low-frequency transport. The overlap between J and 
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the slow modes is captured by the off-diagonal elements of the conductivity matrix. In practice, only a small number 
of conserved quantities is included in the set of slow modes, but the approach can be systematically improved since 
adding more conserved quantities increases the lower bound for the conductivity. 36 

It is convenient to invert the Kubo formula for the conductivity matrix using the projection operator method. 36 
We introduce the scalar product between two operators A and B in the Liouville space 

(A\B) = y [ 1/T dr(A^e Hr Be- Hr ). (F7) 
L Jo 

The susceptibility matrix x is defined as 

Xnm=T- 1 (J n \J m ). (F8) 

We denote by 

V = l-T- 1 J2XnL\Jn)(J m \ (F9) 

nm 

the projector out of the subspace of slow modes. Using identities for the projection operator, equation (F6) can be 
brought into the form 2 

(Turrit) = i{[uO - Mx _1 ] _1 X}nm, (F10) 

where M is the memory matrix given by 

M nm = T~\j n \CV ]— PC\J m ). (Fll) 

CO — / JL, I 

Here, C is the Liouville superoperator defined by CJ n = [H,J n ]. For simplicity, we assumed that all the slow modes 
have the same signature under time-reversal symmetry. It can be shown that if there is an exact conservation law 
(local or nonlocal) involving one of the slow modes, the memory matrix has a vanishing eigenvalue, which then implies 
a finite Drude weight. In the following we apply the memory matrix formalism to calculate the conductivity for the 
low-energy effective model (Al) at half-filling, allowing for the existence of a single conserved quantity Q, [Q,H] = 0. 
The conductivity matrix is then two-dimensional. We choose J\ = J and J2 = Q± = Q — J(J\Q)(J\J)~ 1 so that 
X is diagonal. At low temperatures T <C J, we can use the current operator in equation (E6); to first order in A±, we 
obtain 

where b\ is defined in equation (B7). Likewise, X22 = ((Q 2 ) — (Q J) 2 / (J 2 )) / (LT) can be calculated within the low 
energy effective model once a conserved quantity Q has been identified. The remaining approximation is in the 
calculation of the memory matrix to second order in umklapp. This is analogous to the calculation of the self-energy 
n^ et in equation (B4). Using the conservation law, we can write 

MkM u {w)( * (F13) 



where r = (QJ)/{J 2 ) and 



MnH^^^t-^-^), (F14) 



with 62 as given in equation (B7). Thus, from equation (F10), we find 



where 



y i 1 
+ 



1 + yoj 1 + y to - (I + y)xn M u (lu) 



(F15) 



y=^ = W£ . (F16) 

y X22 (J 2 )(Q 2 )-(JQ) 2 { ' 
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Equating equations (F4) and (F15), we find that the memory matrix approach is equivalent to adopting the self-energy 



loXuM u (lo) -hu 2 = 

1 - yxuMu(u>)/u> ~ 1 - yxIiMni^/w 



-buj 2 — 2ijuj 



(F17) 



As expected, the memory matrix result reduces to the self-energy result for y — > 0. The difference between the 
self-energy and the memory matrix result is of higher order in the umklapp interaction. Therefore, the conservation 
law is not manifested in the lowest-order calculation of the self-energy. Although equation (F17) is not correct beyond 
0(A 2 ,A±), it suggests that the memory matrix approach corresponds to a partial resummation of an infinite family 
of Feynman diagrams which changes the behavior of H(oo) in the limit uj — > from II — > —2i^uj to II — > — y~ 1 uj 2 . 
It follows from equation (F2) that 



lm(J;J) ret (uj) = -ujRea(uj) = -uj 



vK 
2^ 



1 + 2/ 



S(uj) + 



2 7 



(1 + + 6' 2 ) 2 u; 2 + (2 7 ') 2 



(F18) 



where b' 2 



(1 + 2/) 62 and 7' = (1 + y)j. The first term on the right hand side of equation (F18) can be associated 
with the ballistic channel and the second one with the diffusive channel. Substituting equation (F18) into equation 
(Fl) and computing the integral, we find that for times t ^> (2ttT)~ 1 



C(t) 



vKT 
2tt(1 + y) 



-2-f't 



1 + 61 + v 2 



Therefore, in the limit t — > 00 the current-current correlation function approaches the value 

vKTy(l-h) (J 2 )y _ (JQ) 2 



lim C(t) 



2tt(1 + J/) L(l + 2/) L<Q 2 )' 



(F19) 



(F20) 



consistent with the Mazur bound for the Drude weight. 

For intermediate times (2ttT)~ 1 <C t <C 1/7', we obtain the linear decay 



C(t) 



2tt(1 + 6) 



(1 - 2t*), 



(F21) 



independent of ?/ if 61, 6' 2 <C 1. Therefore a small Drude weight cannot be detected in this intermediate time range. 
However, we can extract the decay rate 7 by fitting the numerical results for C(t) at low temperatures (see table in 
fig- 3). 

In addition to the data shown in Fig. 3, we also calculated C(t) at relatively large magnetic fields and various 
temperatures. As shown in Fig. 5 we find that in this case C(t) appears to converge to a finite value within fairly 
short times. Furthermore, we find (see Fig. 6) that at large A and large magnetic fields the Mazur bound (E5) almost 
completely exhausts the Drude weight at high temperatures. This is consistent with the findings in ref. [6] which were 
based on exact diagonalization. 



APPENDIX G: COMPARISON WITH QUANTUM MONTE CARLO 

It follows from equation (B3) that the dynamical conductivity cr(g, uj n ) defined in ref. [5,37] is given by 

/ x u n . vK uj n 

(T(q,UJ n ) = a'XretW^n) = fl < h\ 2 < fl < \ 2 2 i O ' ( G1 ) 

<T 27r (1 + b)uj* + (1 + c)v z q z + 2^uj n 

Here, uj n are Matsubara frequencies. The function <j(q,uj n ) can be calculated numerically by quantum Monte Carlo 
(QMC). By fitting the numerical results, one tries to extract a possible nonzero relaxation rate 7. If 7 = 0, one 
concludes that the Drude weight is finite. Fig. 7 compares the numerical results for a(q,uj n ) obtained by QMC in 
ref. [5] to the field theory result in equation (Gl). We also show the result that would be obtained by setting 7 = in 
equation (Gl), as assumed in ref. [3]. Since the agreement is good for both 7 = and 7 7^ taken from equation (B7), 
we conclude that the QMC results are not useful to decide if the relaxation rate vanishes in the integrable model. 
This is expected because this method cannot resolve a relaxation rate which is smaller than the separation between 
Matsubara frequencies, cj n +i — uj n = 2ttT. Even at the Heisenberg point, where 7(T) ~ Tj ln 2 ( J/T), the relaxation 
rate predicted by field theory seems to be too small to be detected by this QMC method. 
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FIG. 5: Current-current correlation function C(t) for A = 0.4 and a magnetization (m) = 0.3. At low temperatures the 
asymptotic value seems to be reached almost instantaneously. At high temperatures the data seem to be consistent with a 
simple exponential decay to a finite value without oscillations. Using an exponential fit and extrapolating to infinite temperature 
we find C(T -> oo,t -> oo) « 0.058. On the other hand, the Mazur bound, eq. (E5), yields 2TL> Mazur = 0.013. 
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FIG. 6: Current-current correlation function C(t) for A = 4.0 and a magnetization (m) = 0.25. In addition, linear fits of the 
data in the regime Jt > 4 are shown. We find that the Mazur bound (E5) exhausts 95% of the Drude weight for T — > oo. 



APPENDIX H: COMPARISON WITH BETHE ANSATZ 

A finite Drude weight at finite temperatures has been obtained in two Bethe ansatz calculations. 4 ' 10 Both results 
use the finite temperature Kohn formula 1 which relates the Drude weight with the curvature of energy levels with 
respect to a twist in the boundary conditions. Zotos 4 calculated the curvature by an extension of the traditional 
Thermodynamic Bethe ansatz (TBA) based on magnons and their bound states (strings) while Benz et al 10 used an 
alternative formulation based on spinons and anti-spinons. The obtained results disagree. Problems with the approach 
used by Zotos have already been pointed out by Benz et al. For the approach based on spinons and anti-spinons, 
the authors point out that "this approach is phenomenological and for the computation of the Drude weight it can 
not be expected to be as successful as for the free energy". Indeed, the authors show that this approach fails at high 
temperatures where exact relations are violated. 

In Fig. 8 we compare the Bethe ansatz results from ref. [10] at low temperatures with the Drude weight we obtain 
by setting 7 = in Eq. (Gl) 
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FIG. 7: QMC data (symbols) for T/J = 0.1 and a system size L = 128 for (a) A = 1/2, and (b) A = l. 38 In comparison, the 
field theory result in equation (Gl) (black solid lines) and the result obtained by setting 7 = (red dashed lines) is shown. The 
agreement is good, in particular, for A = 1/2. In this case it is clearly impossible to decide whether or not the Drude weight 
is finite based on the QMC data. For A = 1 the effect of 7 is largest, however, here also the finite size corrections are large so 
that it seems again practically impossible to detect a finite Drude weight by QMC. 




The agreement is excellent. This Bethe ansatz calculation apparently does get the temperature dependent parameter 
b right while, on the other hand, completely missing the decay rate 7. We want to stress again that our numerical 
results for C(t)/2JT shown in Fig. 3(c) are already smaller than D(T) obtained this way, proving that 7 is nonzero. 
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